NEW  YORK  UNIVERSITY 
INSTITUTE  OF  MATHEMATICAL  SCIENCES 


AFCRC-TN-60-441  LIBRARY 


25  Waverly  Place,  New  YoA  3,  N.Y^ 


^  %     NEW    YORK    UNIVERSITY 


Institute  of  Mathematical  Sciences 
Division  of  Electromagnetic  Research 


RESEARCH    REPORT    No.    CX-48 


A  Hartree  Self  Consistent  Method  for  the  Scattering 
of  Positrons  by  Hydrogen  Atoms 


RICHARD  S.  RUFFINE 


^  Contract    Nos.    AF    19  (6  04)4555 

I    .'^        DA   30-069-0 

^  APRIL,    1960 


t       ;         DA   30-069-ORD-2581  ,    Project    No.    2360 


NEW  YORK  UNIVERSITY 
Institute  of  Mathematical  Sciences 
Division  of  Electromagnetic  Research 
Research  Report  No.  CX-^8 


A  Hartree  Self  Consistent  Method  for  the  Scattering 
of  Positrons  by  Hydrogen  Atoms 

RICHARD  S.  RUFFINE 


V..O^" 


-e^ 


Richard  S.  Ruffi 


Morris  Kline,  Director 


The  research  reported  in  this  paper  has 
heen  jointly  sponsored  by  the  Geophysics 
Research  Directorate  of  the  Air  Force 
Cambridge  Research  Center,  Air  Research 
and  Development  Command,  under  Contract 
No.  AJ'19(604)il555,  and  the  Office  of 
Ordnance  Research  under  Contract  No. 
DA-3O-O69-ORD-258I,  Project  No.  256O. 


April  i960 


il 


Table  of  Contents 

Page 

I  Introduction  1 

II  Formulation  of  the  Problem  6 

2.1  Form  of  the  wave  function  6 

2.2  Variational  method  11 
2.5  The  differential  equations  l6 

III  The  Dlpole  Approximation  21 

3.1  Behavior  near  the  origin  21 

3.2  Asymptotic  behavior  of  y  and  v  23 

3.3  Multipole  potentials  29 

IV  The  Adlabatic  Approximation  33 

V  Results  and  Discussion  39 

5.1  Numerical  results  for  s-vave  scattering  39 

5.2  Discussion  of  results  ^5 
Appendix  I  ^8 
Appendix  II  58 
References  6h 


Ill 


Acknowle  dgement  s 

I  shoTJld  like  to  thank  Professor  Sidney  Borowltz  for  suggest- 
ing this  problem  and  for  his  aid  and  encouragements  during  the 
formulation  and  solution  of  it.   I  should  also  like  to  thank  the 
staff  of  the  AEC  Computing  Center  of  New  Yoit  University  for  the 
aid  they  offered  while  the  numerical  work  was  being  carried  out. 


-  2 


an  effective  potential  for  the  scattered  particle.   Such  an  assumption  is 
probably  valid  at  high  energies  but  then,  the  wave  function  should  contain 
terms  which  allow  for  excitation  and  thus  a  single  term  in  the  Hartree 
series  would  not  be  valid.   At  energies  below  excitation,  the  incident 
particle  is  moving  sufficiently  slowly  compared  with  the  motion  of  the 
boiind  electrons  that  one  would  expect  the  atomic  electrons  to  be  able 
to  follow  the  motion  of  the  scattered  particle  while  it  is  in  the  vicin- 
ity of  the  atom . 

Refinements  on  the  one  product  scheme  have  been  proposed  by  Seaton"- J 
and  Temkin  ^  -^  .      Seaton  suggests  that  the  total  wave  function  be  expanded 
as  a  series  of  terms  consisting  of  solutions  to  the  bound  state  problem 
with  coefficients  which  depend  on  the  coordinates  of  the  scattered  elec- 
tron.  However,  since  one  can  expect  an  infinite  number  of  terms  the 
series  must  be  truncated  in  order  to  compute  the  coefficients.   Having 
once  computed  the  coefficients,  there  is  no  machinery  for  going  back  and 
improving  the  estimate  of  the  bound  state  distortion.   Temkin  uses  the 
adiabatic  approximation  to  compute  the  distortion  of  the  bound  state.   He 
assimies  that  the  incident  electron  is  moving  so  slowly  that  one  can  consi- 
der the  atom  to  be  in  the  field  of  a  static  charge.   By  means  of  pertur- 
bation theory  he  is  able  to  compute  a  distorted  bound  state  wave  function 
which  depends  parametri cally  on  the  position  of  the  incident  particle. 
The  bound  state  function  is  then  used  to  average  the  Hamlltonian  over  the 
atomic  electrons  to  provide  a  one  particle  Hamlltonian  for  the  scattered 
particle.   Again,  this  technique  suffers  from  the  disadvantage  that  it 
does  not  contain  any  machinery  for  improving  the  estimates  of  the  bound 
state  distortion. 
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Our  aim  has  been  to  develop  a  method  whereby  the  self-consistent 
features  of  the  Hartree  approach  can  be  incorporated  into  a  form  for 
the  wave  function  which  will  contain  correlations  and  satisfy  the 
asymptotic  boundary  conditions.   We  have  accomplished  this  by  adding 
a  siMi  of  products  of  one  particle  functions  to  a  product  of  the  unperturbed 
bound  state  wave  function  and  a  function  which  will  have  the  proper  asymp- 
totic form  for  a  scattered  particle.   By  requiring  the  terms  in  the  sum 
to  vanish  at  large  distances,  we  have  been  able  to  obtain  a  set  of  coupled 
differential  equations  for  the  one  particle  functions.   We  can  show  that 
these  differential  equations  lead  to  an  effective  one-particle  Hamlltonian 
for  the  scattered  particle  in  which  the  potential  is  made  up  of  a  svim   of 
induced  multipole  potentials.   We  can  also  show  that  with  a  suitable  set 
of  approximations^  the  equations  reduce  to  the  adiabatic  approximation. 

Although  our  ultimate  goal  is  to  apply  these  equations  to  the  prob- 
lem of  electron  scattering  from  atoms,  we  have  restricted  ourselves  to  the 
scattering  of  a  positron  by  hydrogen.   By  so  doing,  we  need  not  consider 
effects  of  exchange  which  must  add  an  additional  complication  in  the  elec- 
tron problem,  and  we  have  an  exact  bound  state  eigenfunction  which  we  would 
not  have  for  any  other  atom.   Although  there  is  no  experimental  data  for 
positron  scattering  from  hydrogen,  Spruch  and  Rosenberg  ^  -'  have  computed 
a  lower  bound  on  the  zero  energy  scattering  length  which  can  be  used  to 
estimate  the  efficacy  of  the  approximation.   The  problem  has  some  practical 

significance  in  that  the  cross-sections  for  atomic  scattering  are  necessary 

[9] 
to  compute  reaction  rates  for  posltronium  formation  '--'-. 

There  are  two  previous  theoretical  treatments  of  the  low  energy  scat- 
tering of  positrons  by  hydrogen.   One  Is  by  Spruch  and  Rosenberg  L  -^  and 
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one  by  Messey  and  Moussa  ^    -*  .   Spruch  and  Rosenberg  have  shown  that  the 
Kato  L  -^  method  may  be  applied  to  the  positron-hydrogen  problem^  and 
vill  provide  anupper  bound  on  the  scattering  length  at  zero  energy.   They 
examined  various  classes  of  wave  functions  to  find  one  which  gave  the 
lovest  bound  and  they  found  that  the  effects  of  correlation  were  such  as 
to  make  the  atom  on  the  whole  attractive  to  the  positron  whereas  the 
static  model  would  predict  a  repulsive  interaction.   The  application  of 
this  method  to  systems  which  have  a  bound  state  has  not  been  accomplished 
at  non-zero  energies  and  so  is  of  limited  value  for  the  electron-atom  scat- 
tering problem.  Massey  and  Moussa  introduced  the  correlation  in  a  pheno- 
menological  way.   Thus,  if  one  assumes  the  atom  to  be  a  polarizable  mediijm, 

it  will  acquire  an  induced  dipole  moment  in  the  field  of  a  point  charge. 

o 
This  dipole  moment  is  proportional  to  the  field  which  varies  as  l/r  .   The 

potential  of  a  point  charge  in  the  field  of  a  dipole  varies  as  l/r  times 
the  dipole  moment.   Thus,  the  induced  dipole  potential  varies  as  l/r  . 
Massey  and  Moussa  construct  a  one  particle  Hamlltonian  containing  a  poten- 
tial with  this  asymptotic  behavior  and  choose  a  trial  wave  function  for  use 
with  a  variational  principle.   One  may  object  to  their  work  on  two  grounds. 
The  first  is  that  there  is  no  consistent  argument  leading  to  the  introduc- 
tion of  such  a  potential  and  the  second  is  that  the  trial  wave  function 

they  choose  does  not  have  the  form  which  a  solution  to  the  polarization  po- 

[12] 
tential  problem  should  have  l  -J  .   Th^y  find  that  the  atom,  while  more 

attractive  than  would  be  supposed  from  the  static  approximation,  ia  on  the 

whole  repulsive. 

In  Section  II  we  construct  the  approximate  wave  function  and  use  the 

Kohn  variational  principle  to  obtain  a  set  of  differential  equations  for 


positron  scattering.   In  Section  III  we  shov  that  the  approximation 
resTilts  in  an  induced  miiltipole  expansion  and  we  make  the  approximation 
of  reducing  the  series  to  the  dipole  term.   The  relationship  between 
this  approximation  and  the  adiabatlc  approximation  is  demonstrated  in 
Section  IV.   The  results  of  a  numerical  calculation  of  s-wave  scattering 
are  presented  in  Section  V. 
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II  Formulation  of  the  Problem 


2.1  Form  of  the  vave  function. 

We  make  the  usual  assumption  that  the  proton  has  an  infinite  mass  and 

write  the  Hsmlltonian  in  Hartree  units  vith  'fi  =  ra    =>    m     =e=l, 

e    p 

(2-1)       E  =  K^  +  K^  -  ~-  +   -^  ^ 


r    r   ' 
1  2    12 

vhere 

K,  =  -  I  V^^ 

and 

— ► 

r  =  electron  coordinate 

r  s  positron  coordinate 

The  effects  due  to  the  electron  spins  are  neglected.   Since  the  particles 
are  not  identical,  there  are  no  complications  due  to  the  Pauli  Principle. 

We  shall  confine  ourselves  to  energies  below  the  first  excitation  poten- 
tial of  hydrogen  so  the  boundary  condition  is 

ikr. 


(2-2)       'y{\,%)   ^^100^^1^ 


^2 


2 


wave 


where  ^     n^^i'^   ^^   ^^^  ground  state  hydrogenic  function  and  k  is  the 
number  of  the  positron. 

The  Incident  particle  may  be  assumed  to  be  in  a  superposition  of  orbital 
angular  momentum  states  characterized  by  the  quantum  numbers  I   ,m  ,      Similarly, 
the  bound  particle  may  be  assumed  to  be  in  a  superposition  of  hydrogenic 
states  characterized  by  the  quantum  numbers  n,  ^  ,m_^.   Since  the  Individual  or- 
bital angular  momenta  are  not  constants  of  this  motion,  it  is  more  convenient 


to  use  a  representation  in  vMch  the  total  wave  f^mctlon^  Y(r  ,r  ),  Is 
given  by  a  superposition  of  eigenrunctions  of  the  total  angular  momen- 
tum, ^(L,M|r  -r  ).   Thus, 

(2-3)       nr^^r^)  =  g  A^^^  ^(L^Mlr^,;^) 

■where  the   A  are   chosen  to  satisfy  the  ■boundary  condition,    (2-2). 

Since  the  Interaction,    l/r     ,    does  not  mix  states  of  different  L,M,    the 
Y(L,M|r  , r   )   maybe   computed  independently. 

The  'i'(L,M|r  ,  r   )    can  be  expanded  eis   a  series  of  hydrogenlc  states, 
ijf    .         (i'-,)>    vhlch  satisfy  the   equations 

'^nt^m^^^l^^^V^^-l^^m^^V 


Y.    (fi)   s  spherical  hannonic. 


Thus, 

/    ,        u.     .    1  ij.M:m_  .in_  iw    _  i  r    i    —      x 

"2. 


t^,^2,n     m^,m2  1  2  IT.  2 


xG{L,t^,l^,n\r^)Y^  ^  (fi^) 


\^ 


*— — '         r 


^u       (r^)G(L,e^,Vn|r2)3^,   ,    (n^g)    . 


l^>t^,n       12     ^^1^  -"'-  "'  ^"r2 
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The  L  J  c  ,  (L,M;nL,m_)  are  the  Clebsch-Gordon  coefficients^  the 
radial  functions  G{L,1   ,1   ,n\T   )    are  defined  by  (2-h),    and 


M 
The  y  ^  ,    ,      are  elgenfunctions  of  the  total  aagiilar  momentum.   The  eigen- 

12 
value,  M,  describes  the  orientation  of  the  system  in  space  and  since  the 

scattering  is  independent  of  this  orientation,  G(L,  t,  ^  ip;,n|  r^)  is  inde- 
pendent of  M. 

We  ai^  concerned  with  lov  energy  scattering  by  a  hydrogen  atom  in 
its  ground  state.  The  energy  eigenvalues  of  hydrogen  in  Hartree  units 
are 

(2-6)       ^n  =  "  "^  '    ^  "   1.2,...  . 

2n 

If  the  atom  is  Initially  in  state,  E  ,  the  total  energy,  E,  is 

k2 
(2-7)       E  =  ^  +  E^ 

.    [8] 
where  k  is  the  wave  number.      We   shall  take  n  =  1  and,    following  Spruch    •-    -" 

restrict  ourselves  to 

,  2 

k     <.5 

so  that  excitation  to  states  with  n  >  1  and  pickup  of  the  electron  by 
the  positron  are  excluded. 

Thus,  taking  the  limit  of  f (L,M|r  ,r  )  for  large  r^,  we  see  that 
all  the  G's  must  vanish  except  the  coefficient  of  u   (r^),  G(L,0, i2,l|r2) . 
This  coefficient  has  the  asymptotic  form 
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Lrt« 


L3t> 


(2-8)  G(L,0,L,n|r2)  ->    a^  sinCkrg  "  -2")   +   ^l  '^°^(^2  "  T^" 


Note  that   t     =  L  for  ^     =  0.      The  ratio   Pr/c^  is  real  and  is   given  by 


(2-9) 


—  =  tan   ri_ 


vhere   t]     is  the  phase  shift. 

Let  us  nov  split  off  the  u     (r   )    term  from  ^(L,M|r  ^r^)    in   {2-k) 
and  write 


(2-10)  ni^,U\\,r^)   =  ^  u^Q(r^)G(L.0,LAlr2)y  LoL 


12 


1        V 


V2  e^,^^ 


nitl  1 


M 

y  Lt,^2 


=  yL,M|r     r   )   +  ^       Y]       W  ^        (^^^2^^  U   e    * 
^  ^     ^  ''l''2     t^^tg  ^r2     ^     ^  ^  Lt^tg 


Here,  ^  is  the  form  of  the  vave  function  used  in  the  static  approximation 
n 

and 


(2-11)      W^^.Jr^^rg)^  g  u^^^{r^)G(L,l^,l^,n\v^) 
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One   can  obtain  an  approximate  solution  to  the  scattering  problem  by- 
taking  a  finite  number  of  terms  In  the  series   defining  W  and  reqiilrlng 
that  the  'P(L,M|r^  ^r^)   so  defined  satisfy  the  Hartree  equations 

(2-12)         /    A  -n^^(-l)}^U^^,[:^-^]"^^'^l^r^2)^V"2  =  °- 

An  alternative  approximation  which  ve  propose  is  to  replace  the  sum 
over  n  in  W  .  ,  by  a  single  term  of  the  form 

(2-15)    W  \  ^   (^i^i-g)  =  ^i   (^l^i)^^  (^1^2^ 


where  the  v.    (L|r   )   and  the  y.    (LJr   )    are  to  be  determined  by  a  variational 

calcrulation.      If  we  think  of  W        ,      as  a  correction  to  the  lowest  or- 

\^2 

der  Hartree  term,    f   ,    then   (2-13)  niay  be  thought  of  as  a  Hartree   approxlma- 

H 

tion  for  the  correction.   This  approach  suffers  from  the  usual  difficulty 
in  that  the  correlations  we  obtain  by  summing  over  n  are  suppressed.   On  the 
other  hand,  the  total  wave  function  for  given  L,M  is  coirelated. 
We  define  our  approximate  wave  function,  f   ,    as 

(a-lll)    »p(L,M|?j^,r2)  <.  ■^Y^^{r^)G{h,0,-L,l\r^) 


Z   v^jL|-x)7j^(Llr,)3/»j^,^ 


Note  that  the  stm  over  J  and  J  is  restricted  only  by  the  condition 
J^  +  Jg  g  L  g  1j^  -  Jgl- 
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The  boundary  conditions  on  the  y     (L|r   )   are 


y     (L|r   )  — >0 
^2  "^  r^ 

T  (L|r  )  must  vanish  at  the  origin  and  be  regular  over  elLI  space. 

The  choice  of  functional  form,  for  v  (L|r  )  and  v  (L|r  )  Is  still 

riot"^  2  r  -| 

an  open  question.  Massey  and  Moussa  ■-  -■   and  Spruch  and  Rosenberg  L  ^ 

Introduced  analytic  expressions  for 

These  expressions  contained  unknown  constants  which  were  determined  by 
a  variational  principle.   Such  a  procedure  enabled  them  to  include  terms 
with  an  r   dependence.   Having  defined  W  .  .  (r  , r  )  to  be  a  product 

-Lf^  J  -I  J  p        -L        ^ 

of  one-particle   functions,   we   are  free  to  determine  these  functions  by 
means  of  a  variational  principle. 


2.2     Variational  method. 

In  order  to  obtain  differential  equations   for  G,   y.     and  v       we 

^  Vll 

shall  use  a  variational  principle  of  the  type  Introduced  by  Kohn    L  J  , 

We  first  simplify  the  notation.      We  shall  suppress  the  explicit  radial 

Dependence  of  G,v  and  y  as  well  as  the  explicit  L  dependence  of  v  and  y. 

Then^   we  may  write 
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G(L,0,L,l|r2)  =    G^ 


^J.^^1^2)  -    ^J, 


We  Biay  also  'vrrlte 


^lo^^i^  =  "  °  ^V 


With  this  notation,    (2-lii-)   becomes 


(2-15)  ^L^Mlr^^r^)   =  ^p  =  T^T 


L  JLoL  "^     -f^ 
''1*^2 


uG.y 


.M 


^j/j^y^JiJ, 


The  first  order  variation  in  f     is  defined  hy 


(2-16)         5>P    =  — r 


"^•^/loL 


and,    afl  r  — >  » 


5y       ->  0 


5v       ->  0 
"^1 


(2-17)         6Gj^    ->  5a^  8in(kr  -  ~)  +  ^^^  cos(kr  -  -g-) 


The 


first  order  variations  given  in   (2-1?)    insure  that  the   function, 


f     +  5¥  ,   vill  satisfy  the  "boundary  conditions. 
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We  define  the  Integral, 


(2-lQ)    I  =  /  '^p  (H-E)>?pdT^dT2, 


and  compute 


(2-19)    51  =  r5'Fp*(H-E)7pdT^dT2  +  J  ^^*(E-E)bY^i-:^d-x^ 


The  remainder;,  A  ,  arises  farom  the  integrations  by  parts  vhich  are 
necessary  to  bring 

r¥p*(H-E)5¥pdT^dT2  into  f  bl^{E-E)l^\x^d':^  . 

When  these  involve  "u,v.  ,y.  ,5v   and  5y   no  remainder  accrues  ovlng 

J-L  ^2   "^1       "^2 

to  the  vanishing  of  the  boundaiy  terms.   The  integrations  by  parts 

necessary  to  transpose  G,.   and  5G^  are  carried  out  by  taking  the  limit 

■Li         L 

of  the  integral 

R  sr  * 

1^*1       d         2d         L        2^  TllS, 

r^     ^; fe   r^     dr„      /      u       dr,     X 


(2-20)    -i/ija 


^     L         2  dr„     2     dr_       2      2         2     /     '    '  1 

r  r  2  2    r  ^ 

o  d  o 


><    /^LoL>^!oL^^^2 


Ik 


R 


1   r    *  d 

-'      dr 


as  R  — 5»-  oo.   Integrating  by  pariB  tvlce  leads  to 


R 


(2-21)   -i 


*  d 

^L  -2  5^ 
dr 


*  d 
G^  ~   SG  - 
L  dr   L 


d    * 


r=R 


R      p 

1   r       d     * 


The  last  term  of  (2-21)  supplies  the  desired  transposition.  The 
bracketed  term  may  be  evaluated  In  terms  of  the  asymptotic  forms.   Thus, 
using  (2-8)  and  (2-17),  we  have 


(2-22)    A^=-| 


'  *  d      ^\ 

0,  I-  6G  -  -—-   &G^ 

Tj  dr   L   dr    L 


R 


=  -iKPLX-^X 


Recalling  that 


°1 


;^)=tan,^. 


A  may  be  written  as 


(2-25) 


^1=  -  2 


tan  r|^  5(a^  a^)    -   5(a^  ^) 
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We  assume  that  5v  ,    5y   and  6G  are  Independent  variations  and 
ve  further  assume  tan  ^  in  (2-23)  to  "be  set  at  its  correct  value.   Then, 

(2-21+)      51  -  A^  =  5  Tl  +  |(a^\  tan  n^  -  o^^P^)"^  =  0 


if  the  coefficients  of  the  independent  variations  Sv  ,  5y  ,    5G  and 
their  complex  conjugates  vanish.   This  leads  to  the  stationary  expression 
for  the  tangent  of  the  phase  shift;, 


(2-25)      I  -  I  o^^^   =  -  I  a^\  tan  tIj 


Requiring  the   coefficients  to  vanish  leads   to  the  equations  defining 
V      ,    y        and  G^^;    i.e. 


(2-26)  /     ^^yl^^   (H-E)^pdT^dn2  =  0 


(2-28)  Z    /    ^  -jj  y  LoL^«-^)V^1^2  =   ° 


Jl 


and  a  similar  set  for  the  complex  conjugates. 
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2.5     The  differential  equations 

On  performing  the   Indicated  Integrations  In  (2-26^27,28)  we  have 


(2-29) 


K^^)   +  V. 


H  2- 


G^-H    <ulv^> 


'M.     1        k^ 
^  '  "2 


Il<-\\\^,^<yl^\\\¥\,^^-,, 


J-j_J2 


(2-50)  N 


(Jl) 


(Jl)       1 
K  -  — 


(Jl) 


r,  2 


z 


<^J  l''rl>'3i> 


2''       2''    I'' 


<J'LJ,J,l^l^LJiJ.>^Ji*<^L 


°L^  "'j.c 


^<\\\<^'><f\>,\\flL>-'° 


(2-51)  N 


(Jg) 


(Jg).    1 


(Jo)   1 


K  + 


^ 


<      -I    •^2     7,J^,J{,J^  -^l  '^l 


<yLj,jj^7iyLj.j.>-j.^<^oi- 


> 


^2         2 


^L  ^J,L 


V^ 


M  1^      i^,M 


2.    <'jJ«jJ"><KvJ^jJKoL>«L 


=  0 


-  17  - 


vhere 


(2-52)    K^JL-l^.lLi!^ 


2  ^  2 
dr 


2r 


We  may  ignore  the  complex  character  of  the  radial  functions  since  it 
enters  only  through  a  complex  normalization  on  G  and  consequently  may 
be  cancelled  out  of  the  equations.   Thus, 


(2-55)     <v|f|z>  =  r  w*(x)f(x,y)z(x)dx 


=   /  w(x)f(x,y)z(x)dx 


The  integrations  over  l/r  p  ^^'^  carried  out  by  expressing  l/r   as 


(2-3^) 


"12     7 


< 


y       (r>) 


7+1  ^7^^°^  ^12^ 


=  y 


R^P^(co8  9^2 ) 


-<=  < 


r^>    Tg  >  r^ 


r^,    ^2  <  ^1 


vhich  defines  R, . 


^' 
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The  angular  integrations  are  represented  by 


^^-55)     Q\^^K\f\y^y    -    /^LJ^j/7^-«  ^12)^LJ'J'^l^V 


[6J 
These  coefficients  have  been  tabulated  by  Seaton     for  a  number  of 

values  of  J  and  J  .   They  are  Independent  of  M  and  vanish  unless 
(2-56)    (-)  1   2=  (-)  ^ 


(2-37)    (-)  ^   ^    =  1 


(2-58)    j^  +  ^  ^  ^2  -  1^2"  ^1- 

Equations  (2-36,57)  are  expressions  of  parity  conservation  and  (2-58) 
means  that  states  vith  different  multipole  symmetry  cannot  be  mixed  by 

a  multipole  of  order  lower  than  the  difference  betveen  their  orders  or 

M 
by  an  order  higher  than  their  sum.   The  V  y  ^  j,  satisfy  the  orthogonality 

relation 


V_  is   defined  as 


v,.i-  T.  <A^\-><vl,j.K\y-t^> 
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vhlch  by  (2-38)    and   (2-39)    reduces  to 


(2-4o)        \-^-  <^\-rJ^> 


1 

r 


i  r.^a. 


2 


e-^^(i-M). 


The   constants  E     and  E     are  defined  by 
V  y  "^ 


(J.)  v-n 

(2-1+1)  E^     =  -   2     2^     <: 


J2  2 


(Jg)      1 

K  +  i  -  E 

r 


'i,""/  ^ 


(Jl) 


(Jo) 


(2-1+2)  E^     =   -   2     2-j      <■ 


h  "^ 


(V    1 

K  +  -  -  E 

r 


/(J2) 
V,    >/n 
J2     /      y 


E=|(k2  -   1), 


(Jl) 


(2-1+3)       N^"-   =   2-   <yj 


|y.  > 

J2  "2     ^2 


(Jp)  v-^ 


^  ^1      ^     "^^      ^^ 


The  B-tans  over  J^  in   (2-1+1,1+3)    are  restricted  "by 


L  +  J^  g   Jg  ^    |L  -   J^|. 
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A  slallar  restriction  holds  for  the  sum  over  J  in  {2-k2,kk) . 

In  general,  E  and  E  will  not  be  eigenvalues  of  (2-50,51)  and. 
consequently  the  complementary  functions  of  these  equations  vill  not 
satisfy  the  boundary  conditions.   Thus,  v   and  y   must  be  the  parti- 
cular  integrals  of  (2-50,51).   There  is  the  possibility  that  for  suffi- 
ciently large  k,  E  will  be  positive  and  (2-51)  will  have  a  complementary 
function  with  the  asymptotic  form  Bin( >/e~  r) .   Such  a  term  must  not  be 
thought  of  as  a  possible  way  of  introducing  excitations  since  in  general 
the  solution  to  (2-50),  parti culaj"  integral  or  otherwise,  will  not  be 
aji  eigenfunction  of  the  hydrogen  atom  Hamiltonian. 
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III     Ttie  Dipole  Approximation 

5.1     Behavior  near  the   origin 

The  differential  equations   (2-29_,  50, 51)    and  the  boundary  conditions. 


(3-1) 


y,    ->o 
J2 

,      r  ->  0,00 

V       ->  0 

■^1 

,      r  ->  0,00 

^L     -^^ 

sln(kr  -  ^ 

^L     ^'' 

r  ->  0, 

l^  cos(kr  -  —) 


define  an  infinite  set  of  v,  's  and  y  's  for  each  value  of  L.   These 
equations  must  be  reduced  in  number  before  a  practical  solution  is  pos- 
sible. We  shall  accomplish  this  reduction  by  assuming  that  the  most 
Important  contribution  of 

J-iJp       W-]   up 

to  the  phase  shift  is  made  in  the  region  of  large  t^.     We  shall  then 
shov  that  the  leading  term  In 

and  that  this  term  gives  rise  to  an  effective  Induced  dipole  potential 
in  the  equation  defining  G  ,  (2-29).   This  approach  may  be  Justified 
by  observing  that  in  the  neighborhood  of  the  origin  the  behavior  of  G_ 
is  determined  primarily  by  the  homogeneous  part  of  (2-29).   To  see 
this,  we  note  that  (Appendix  l) 
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(5-5) 


<v|R^|z> 


r  <v\-y:^\^>  >  r 


We  use   (5-5)   to  write  equations   (2-29^31)   in  the  limit  of  small  r. 
Thus, 

2  -I 


(5-M 


K^^)+V^-^' 


^L  ^   <"l^o> 


^(L)    .1       k 
K         +  7  -  "2- 


<^l    t    IV^L^  ° 


(5-5)  N 


(L) 


r 


,(L) 


^L  - 


2     <v,J  ^  !->  y. 


+    <vju> 


r        2 


^-    <^ol    t    I^>Gl=°' 


Since  V„  — >  —  -  1  the  behavior  of  G^  near  the  origin  is  given  "by 
H    r  L 


(3  6) 


r    2 


N 


N-  <v  u>   -> 


(^=Ey^ 


where  E     is   a  constant.      A  similar  equation  for  y     can  be  de rived. 
Thus,    the  forms   of  G     aad  y     in  the  neighborhood  of  the  origin  are 

Li  XA 

determined  by  the  singularity  in  K    or  by  the  Coulomb  repulsion  if 
L  =  0.   In  either  case. 


(5-7) 


3t  -^  Xj  C  r^ 
L    ■^— J   n 


n+L+1) 
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The   coefficients,    c     and  c   ,    are  Independent  of  y     or  v  ,    c     vlll  depend 

on  y     and  c   ,    n  >  2     -vrlll  depend  on  y  ,    J   ^  L.      Since  the   centrifugal 
L  n  J 

potential;,  [L(lrH)J/r  ,  euad  the  Coulonib  potential  are  repulsive,  we  may 
expect  that  the  positron  will  not  penetrate  the  atom  too  deeply  and  the 
correlation  term, 

will  not  be  important  for  small  v^. 

5.2     Asymptotic  behavior  of  y  and  v 

When  r^;>rp  — >  m  equations   (2-50,51)   take  the   form 


(5-8)  E   [/(J)   5^j   ^  V^j  ]v^  .  H^ 


where,    for  w,  =  y.. 


(5-9)  =^^^).N^^)[K^^)    -E^^^)] 


(5-10)  V^j=     2^      <-jJ\l-J.><^LJ^iI^llyLjp> 


(5-11)         V^j=Vj^     ;     V.j=0,    i^J±   1 


(5-12)         H^=   <-J«J-><^Lt    l^tli^L^^ 


(5-15)  t  =    |l  -  l|. 
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When  w.  =  v^   we  have 


(5-lM  /^i)=N^i)[:K(^)    -E^^)] 


(5-15)       ^ij  =E  <yj  l\|yj.><^^,j  |pjj^^jj.> 


J2J2       2 


^5-16)         V^j=V^^   ;   V^j=0,    1^J±1 


(5-17)  H^=     4.<yj    l^lI^L>^^LlJ    l^ll^LoL^- 


Jg  '^2     ^      ^  ^   ^2 


Equations   (5-9)   thru  (3-I8)    are  obtained  by  taking  the  leading 

terms   in  each  of  the  expressions  making  up   (2-30;,5l)    and  only  the  lowest 

order  terms  have  been  kept  In  each  element  of  the  mB,trlx,    X^      6     +  V.  .. 

The   case  of  w     =  y^    requires  special  treatment   and  will  be   considered 

after  the  other  cases. 

We  assume  w     to  have  the   asymptotic  fona, 
J 


(5-18)  Wj  -^fj(r)   e^^ 


vhere   f  .(r)    is  an  aBynrptotlc  series  of  the  form, 
J 


(5-19)  f.  =     Xic^^'^^   r  ^J  ""   , 

'^  n 
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Then 


(3-20) 


T.     N^^) 


K 


(i)        E 


(i)-l 


'l/j 


.(1) 


f       (x2.E(i))v^=A(^) 


where  we  drop  all  terms  which  vanish  faster  than  re      .     Using  (5-20)^ 
we   rewrite    (5-8)    as 


(5-21)  X!     AW     =     2    Ea^^ 


i  ^^    ^  J 


'ij   ^  ^iJ  >J   =  ^1    • 


The  matrix.    A,    is  real  and  symmetric   and  may  therefore  be   diagonallzed 
by  a  similarity  transformation.      If  we   consider  V       to  be   a  first  order 
term,    then  A  is  diagonallzed  to  first  order  by  the  transformation. 


(5-22)  X! 


\/jk\^  =  '^'^'Ki 


where 


IJ 


(5-^5)  ^IJ  =  ^IJ   -   ,(J)    .  ,(i) 


The  inverse  matrix,  A  ,  is  obtained  from 


(5-2M    A-J=  E  B^^-^ 


-J  "  fet  "i^  a"^  "^^""^J 


li 


^  °ij  -  TD^IT  • 
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The  solution  to  (5-21)  is  given  by 


(5-25)    V,.  E  A-JAj^w^.  S  A-Jh^ 


1   H     ^    /iJ^J. 


The  asymptotic  behavior  of  v   is  given  by  substituting  (5-17) 

■^1 
into  (5-25)  and  setting  i  =  j^  v   =  v   .   Then,  we  must  have  X  =  -  1 

in  (5-19)  and 

(5-26)    V   -^  f  e"^  . 
•^1      '^l 

The  form  of  the  series  is  not  of  importance.   The  case  of  J  =  1  is  cal- 
culated in  Appendix  I. 

Since  v.  vanishes  exponentially, 


^1     ^1   r  -^    r   -J    ^1   ''l      r    ^1    ^1 

o 


2 

Thus,  V  J,  as  defined  by  (5-10),  vanishes  as  l/r  .   From  (5-12)  we  have 

H.  =  rH.  ^     or  H.  =  —  H.  ^  .   In  either  case, 
i     1+1        i   r  i+1  ' 

(5-27)    V^jHj^H^/r  or  H^/r^ 

so  that  the  second  term  on  the  right  of  (5-25)  vanishes  more  rapidly 

than  the  first  and  we  have  for  i  =  J  ,  w.  =  Y,   , 

'—        Jp    Jp 


""""     %^7iT''i- 
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In  order  to  satisfy  (5-28)  we  redefine  w  In  (5-l8).   Thus,  take 


(3-29)    yj  =  fj  (r)sin(kr  -  ^)  +  gj  (r)cos(kr  -  ^) . 


With  (3-29);,   we  nrost  take  X  -   ik,   where  i  =    \/-T  ,    in  previous   formulas. 
Thus 

(,.50)        a(J)  .  -  4^  (e/'=*  -  .^ 

and  from  ( 3-12, 15;.  25,  30 ) 


2  0^    1^  >  M  M 

^5-51)  yL±i— >    ^     j^(L±e),g   (Lit)    _  ^2      ^^Lt(L±t)l^tl^LoL>^L 


^rL±ti  ^ 
"Tjii  ^L  • 

r 


Comparing   (3-51)   with   (5-29),   we  have 

(5-52)  cJ^''^=Ql±,     ,      ^go  =  *«^W±t     '      ^L±t=-^'^^^' 

where 


The  case  of  t  =  0  is  complicated  by  an  additional  term  appearing 
In  H_  and  by  the  asymptotic  cancellation  of  the  G^ /r  term.   ThuB) 
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(3-33)         \=   < 


2-'2'-><^La.l^2l^L0L>\ 


v„  R 


K 


(L)        k' 


2  -, 


Let  us  write   (2-31)    in  the  limit   as   r     ->  «>  for  Jp  =  L.      Thus 


(3-5^)  N 


(L) 


K 


(L)       !y 


y^  +  <vju> 


K 


(L)        k 


2  -r 


^ 


S     -il<-Jt^l->^^LJLl^2lfLoL>^L-° 


J       r' 


In  the  sum  over  J^    the   term  with  J   =  1  vanishes  by  the  parity  rule^ 
(2-56),    and  the   first  non-vanishing  term  occurs  for  J  =  2.      If  L  <  2 
all  terms  in  the  s\m  vanish  becatise  of  the  vector  coupling  rule,    (2-38), 
In  this   case,    (3-5'^-)   reduces  to 


(3-35)  N 


(L) 


K 


(L)        i 


(L)   , 


y^  +  <vju> 


K 


(L)        k^ 


^ 


0. 


From  (2-29), 


(3-36) 


K         -    2 


<u|v^> 


K 


(L)        k^ 
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Substituting  (5-36)  into  (5-35)  results  in  an  equation  of  the  form 


(3-57) 


K^^)  -  e' 


^L  = 


\  yj^l  +  B_  y^_^ 


G  G^ 


so  that 


(5-58)    y^ 


^'   E 


G      L 

~  7 


The  theory  presented  in  this  section  could  not  he  considered 
rigorous  insofar  as  the  treatment  of  the  differential  operator,  K,    is 
concerned.   HoweTer,  a  more  extensive  treatment  involving  the  Green's 
function  for  the  system  is  not  necessary  since  we  need  only  the  leading 
terms  of  the  asymptotic  series.   If  the  values  of  y.  are   inserted  into 

(2-51)  J,  the  dominant  terms  in  the  differential  equation  are  J-ust  the 

p(±) 
ones  chosen  In    cC^    '    and  V.  ..  We  may  therefore  claim  that  the  solutions 

obtained  above  do  not  lead  to  any  inconsistency. 


5-5  Multipole  potentials 

The  asymptotic  form  of  y   given  in  the  previous  section  will  now 
be  used  to  determine  which  values  of  J  and  J  make  the  largest  contri- 
bution to  the  scattering.   Let  us  first  consider  the  case  of  y  .   Now, 


(3-59) 


■  (L)   k^^ 
K    -  -2 


^L^ 


K^^).  e'+  (e'-  %-)  ^ 


Jr 


\ 


+  (e 


C  G. 


(%-  E^r"" 


^     ^^'7, 


& 


r5 
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so  that  Tising  y     in   (2-29)    contributes   a  term  vMch  vanishes   as  l/r^ 
Ij 

for  the   case   of  L  <  2.      All  other  y,      contribute  terms  vhlch  vanish 
asl/[r2(^-^l)j.  ' 

Substituting   (5-3l)   into  the   last  term  of  (2-29),   ve  obtain 


(5-^0)  E    J2     <-|Hjl-J><5^LoLl^jI^LjlL±^^^,L±t| 


±1  T 


sr     ^i 


4^       2(1+1)    ^L  =  \  ^• 


The  "potential" 


"^M  "     -f-r^       2(^+1) 


is   a  sum  of  Induced  multlpole  potentials.      The  leading  term  Is  given  by 
■1   =  1  and 

r 
The   induced  quadrupole   term;,    1  =  2  supplies  the   first  correction  and  it 
vanishes  as   l/r   .      Thus  we  drop  all  t  >  1. 

When  L  g   2,    y     can  contribute  a  term  to   (2-29)   which  vanishes   as 
G  /r     which  arises   from  the  term  dropped  from  (3-5^)'      Or   vould  therefore 
be  proportional  to 


(5-^1)  <  Vpjt^lu  >     ,      L  g   2. 
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J+1 
It  Is  shown  in  the   appendix  that  v .  — >  r         in  the  neighborhood  of  the 

origin.      Thiis,    as   J   increases,    the  first  maxinnxm  of  v     moves   out  from 

J 

the  origin  and  consequently  the  overlap  of  v     and  u  can  be  expected  to 
decrease  vith  increasing  j.      Then, < v    |  f^  |u  >  vlll  decrease  vith  in- 
creasing  j   and  since  Q,^.    is  proportional  to   <v   |t    |u  >  we   can  expect 


^±1     >     ^t      '      '^^' 


If  these  estimates   are   correct,   we   find  that  the  largest   contribution  of 
W    .    .      to  the  effective  potential  seen  by  the  positron  arises   from 
•£•  =  1,    J     =  L  ±  1.     MaMng  this   choice  for  Jp  restricts  us  to  a  single 
value  for  J^,    J     =  1. 

The   asymptotic  form  of  f     in  this  approximation  is 


(5-^2)  'P^  — > 


P  r^r^ 


M  1  /  M  M  \ 

"  ^LoL  +  7  "l  {^1  Kl  L+1  ^  \-l  fm  L-l) 


^L- 


This   result  is  in  agreement  with  perturbation  calculations  performed  by 
Bransden,    Dalgamo,  John    and  Seaton    '-    -^    and  by  Temkln    ^  -' . 

By  retaining  only  the  terms   for  which  j     =1,    J=L±lwe  effec- 
tively assume  that  the  most  Important   distortion  of  the   atom  arises   from 
the   coupling  Introduced  by  the  dlpole  term.      By  dropping  y     and  v     we 
imply  that  the  spherical  distortion  of  the   atom  which  arises  from  the 
monopole  term  in  the  expansion  of  IJt  ^  is  less   Important.      It  should 
be  noted  that  our  argument  for  the  unimportance  of  y     hinged  on  the 
fact  that  the  leading  term  of  the  Inhomogenelty   <v    |r    |u>G     in   (2-51) 

O        O  J_j 

was   cancelled  at  large   distances  by       -~F:<y    |u>G^.      If  this  had  not 

°  ■'  2       o '  L 

r 
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been  the  case,  y  would  have  had  the  asymptotic  form 

Li 

and  the  Inhomogeneity  In  the  equation  defining  G  vould  have  been  of 

o 

the   form  G^ /r   .      Of  course^    If  the   cancellation  did  not  taice  place, 
the  atom  vould  have   a  net   charge   and  the  effects  of  the  spherical  corre- 
lation might  be   svemped  by  the   Coulombic  term. 

The  arguments  vhich  vere  presented  here  hinged  largely  on  the   form 
of  V.      at  sniall  distances   and  on  the   form  of  y.      at  large  distances.      If 
these   arguments   are  valid,    they  supply  some   justification  to  our  choice 

of  W  ■^    .(r,r-)    in  2-15).      Thus,    the  exact  W   ,    „      is   given  by 
\^2  ^  ^  12 


(5-^5)  W^^     =     X]     ^^,  G(L,^^,^2.n|r2). 


1  2 


In  general,    u   .     behaves  near  the  origin  as 
""  1 

V^        -^   Vi^     • 


So  that  for  small  r  , 


W  \  ,       >     r^  1         2     e^G(L,t^,l2.n|r2) 

1  2     r^   — >  o  n 

and  If  this   is  the  only  region  of  importance  vhere  the  form  of  W   ,    , 

V2 
is   concerned,   ve   are  then  Justified  in  assuming 


\^2         ^1     ^     ^2     ^ 
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IV     The  Adiabatlc  Approximation 


:n  W 


Temkin    ^  -I    obtains   the   adlabatic   approximation  by  assuming  that 
the   atom  is  in  a  static  field  whose  source  is  the  incident  particle. 
The  static   field  is  treated  as   a  perturbation  and  a  distorted  bound 
state  function  is   found  which  depends  parametrlcally  on  the  position 
of  the  incident  particle.      The  total  wave   function  is   a  product  of 
this  distorted  bound  Soate   and  a  function  which  depends  on  the   coor- 
dinates of  the   scattered  particle.      The  total  wave   function  as   com- 
puted by  Temkin  for  electron  scattering  from  hydrogen  is 


(^-1) 


=  ^^^1^   ^AD^^l^    P^(cos0^2h    ^("^^l) 


loo 


(^«) 
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where 


e{T^,T^)   = 


0   ,    1-2  <  ^1 


1   .    ^2  -^  ^1 


P  ( cos  0      ) =    Legendre  Polynomial 


-AD  =    ^2r'  ^   -')^" 


Now,    (^-1)   is   strongly  reminiscent  of  (^-k2) .      Thus 


(5-^2) 


Y. 


00  1     ^ 


-^ 


M 
LoL 


M 


V^^^^giV-^.. 


^ 


v(r^) 


=    V, 


^ 


Equation   (^-l)    can  be   obtained  from  (3-^2)   by  assuming  that  y        is 

nea;lia;lble   for  r„  <  r     vhere  r     is  acme    cut-off  radius   and  that  for  r„   >  r 
^     ^  2  c  c  2  c 


\ 


(i^-2)        y.    =  Q.     2 

^2  ^2  r^ 


vhere 


(^-5)  Q 


<  V    r    u  > 


M        1       I      M 

<y,, ,   Ipj  vr„r  > 


'J2'  ,Th?^ri7X|v>  ^^LiJs'  1'^^°^ 

I      r   21 


We 


must  compute  v(r  )  for  r  <  r  and  to  do  this  we  assume 


{^-^)  <y. 


r+i 


y   >  =0,7^1. 


Applying  {h-k)    to  E  ^  '    results  in 


(1) 


(i^-5) 


2   "2 


since 


K 


(J2)   k^- 


1  dy 

V   <v/  —  — — 


The  differential  equation  for  v     is 


(2-50)  N 


(1) 


K 


(1)         1        \ 


r  2 


Jg./.J^  ^2  ^2 


35  - 


'^  <  S^Llj/ll  J'loL^  -  =  °- 


By  (^-^),  we  may  drop  the  second  term  In  (2-50)  since  ve  are  interested  in 
r  <  r  .   The  last  term  in  (2-50)  is  reduced  by  (^-2).   Thus 


(^-6)        <yjjHjG^> 


=  -1    /    yj,  7  ^-2 

T  -^2 

c 


J2     r 


^j/j/"2 


Jo 


Substituting   (^-6)    and   (^-5)    into   (2-50)^,    we  have 


(^-t; 


V(l)        1  +  1 


V     =   r  u 


^       1,^(1)         1  1|       V 

<VK^       ' .    +    7rV> 

I  r       2' 


<  V  r  u> 


or 


;^-8) 


d  2  2        , 

:;^-^+ r-  ^ 

■  dr  r 


V  =  hru 


where 


<  V 


(1^-9)  h  = 


r-  2 

"^    d  2  2, 

7^--  +  --  1 

dr  r 


v> 


<v    r    u  > 
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Nov,    h  is   completely  arbitrary  since 


2   .  ,.M 


,M 


<^i---i<n^\\\Yi^^> 


and,  since  the  solution  of  (^-8)  is  directly  proportional  to  h,  the  pro- 
duct, "V-,y.  y    is  independent  of  i 
solution  to  (ij-S)  for  h  =  -2  is 


duct,  V  y  ,  is  independent  of  h.   It  is  shown  in  Appendix  I  that  the 


(1.-10)    V  =  v^. 

The   coefficients  <i^  |P    |  ^         >     are  evaluated  by  Seaton    LJ 


Thus 


/ 


v/L 


(^-11)  <^^',,JP,I^L>     ■     < 


n/5(2L  +  1) 

-  v^L  +  1 
L    V5(2L  +  1) 


,      jg  =  ^   -   1 


,      J2  =  L  +  1 


Let  us   consider  the   case   of  L  =  0.      Then,   with  h  =  -2, 


(^-12)  Q 


'J, 


v/3 


,   L  =  0 


and 


(^-15)  2/°^=       X       C^l^0^0;M^,m2)Y       Y^^ 

m^,  m^  1       <= 


-  57  - 


C^^(0,0;m^,-H^)Y^Y^_^ 


-^  lfeV^°^^12) 


Thus 


P     r,  r 


12  I- 


uU  +  -pr  v„^  T—  P-(cos  0_-) 

Tooo    2  AD  ^Tt  1      12' 


or 


T 


^oo  ^  -2 
^2 


1   ^AD  ^1^"°^  ^12^ 


(^rt)- 


G 


r^{h^) 


1/2 


Comparing  (^-1^)  vitli  (^-l)  we  see  that  the  equations  are  the  same  except 
for  the  difference  in  sign.   This  results  from  the  fact  that  Temkin  con- 
siders electrons  while  we  have  considered  positrons. 

The  case  of  L  >  0  cannot  be  examined  without  a  more  detailed  descrip- 
tion of  the  function  G(r  )  in  (^-l).   It  should  be  noted  that  one  of  the 
features  of  our  method  of  calculation  is  that  angular  effects  are  treated 
systematically  throughout  the  development  of  the  theory.   This  is  not  the 
case  with  the  adiabatic  approximation  as  derived  from  perturbation  theory. 

We  note  that  the  derivation  of  v   is  independent  of  L.   Thus,  we  can 
use  V   and  Q.  to  compute  the  last  term  in  (2-29),  the  equation  defining  G^ 
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We  have 


(2-29) 


K<-)   .  V„   -  I 


2  -1 


G^-     X      <^K|vi>    <yLoLl^lli^Llj? 


Lljg'    ^Jg  =  °" 


Here,    j     =  L  +  1.      Nov, 
(1.-15)         <u|Ri|v^>    =^1 


ut  v^dt  +  r     /      u  -2  v^  dt    . 

r    ^  ^  t 

o  r 


We  note   that  v        Is   regular  over   all  space   and  since  It   decays   exponentially 
ve  may  approximate   (4-15)   by 


(^-16)  <u|Rjv^>_>-^     <u|rjv^>     =2^ 


r  -^00  r 


Using   (4-2),    (4-11)    and   (4-l6),   we  have 


(4-17; 


k(-)  .  V,  -  ^ 


a     L 

'L        2  "T 
r 


a  =   4.5 


for  large  r.   The  last  term  in  (4-17)  is  an  induced  dipole  polarization 
potential.   The  value,  a  =  4.5,  is  a  well  known  result  for  the  polarizlbility. 
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V     Results   and  Discussion 

5.1     Numerical  results   for  s-vave  scattering 

Equations  (2-29^  50_,  31)  were  evaluated  numerically  for  s-vave  scat- 
tering at  energies  such  that  0.  S  k  g  .6.  The  dipole  approximation  vas 
assumed  vhich  reduced  y       to  the  single  function^    y   .      With 

G     =   G 
o 


(5-1)         y^  =  y 


T^   =   V       , 


ve  have 


(5-2) 


dr 


^  -   2Vjj  +  k^l  G  =  -^<ulR^|v>y 


\/5 


(5-3) 


2 

d  ^        2       _2_ 

^2"  2'^r'^N 

dr  r                    V 


<y|Roly>  +  ^<ylR2|y> 


+  E      >v 

V 


V^ 


<y|R^lG^>u 


(5-^) 


d^         2 


^    2  2 

dr  r 


2        2_ 

r  '*'  N 


<v  R     v>    +  -  <v  R     v> 
'    o '  5  "^ 


*^r 


<  V  R-    u  >   V 


^/3^    '^ 


ko 


The  nonaallzatlon  of  Y  'depends  on  the  normalization  chosen  for  G.  For 


k  =  0,  G  was  assumed  to  have  the  form 
(5-5)     G  >  r  +  A^, 


o 


■where  A     is  the  negative  of  the   scattering  length.      For  k  5^  0;,    G  was 
assigned  the   asymptotic  form 

(5-6)  G     >  sin  kr  +  tan  x]   cos  kr. 

r  —>  00 

In  Table  1_,   we  ILst  the  values   of  the  phase  shifts  obtained  by  the 

[81 
self-consistent  method  and  those   found  by  Spruch  and  Rosenberg    '-  -^    ,    and 

Massey  and  Moussa   ^    -^  .      We   also  list  \mder  the  heading  "static  model", 
the  values  that   are  obtained  by  setting  v  =  y  =   0  In   (5-2).      This   data 
is   displayed  in  Figure   1  where  we  have  plotted  tan  -q/k  vs.   k. 

From  Figure  1  we   see  that  the   self  consistent  method  gives   a  positive 
phase  shift  for  k  <  5-5«      At  higher  energies  the  phase   shift  becomes  nega- 
tive  in  agreement  with  the  results   of  the  Born  approximation.      Comparing 
this  with  the   static  model  for  which  the  phase   shifts   are   always  negative, 
we  see  that  at  sarfficiently  small  energies,    the  effect  of  the   correlation 
Is  to  make  the   atom  attractive  to  the  positron.      The  resiilts  of  Massey  and 

rioi 

Moussa  ■-  -'  indicate  that  the  addition  of  an  ad  hoc  polarization  potential 
Into  the  Hamlltonian  will  not  greatly  increase  the  attractiveness  of  the 
potential.   The  values  listed  for  Spruch  and  Rosenberg  '-  -'   are  the  nega- 
tive of  upper  bounds  on  tan  ri/k.   Thus,  In  this  context,  they  represent 
lower  bounds.   The  values  of  tan  r^/k  as  given  by  the  self- consistent  model 
are  less  than  these  bounds  and  must  therefore  be  considered  to  be  in  error. 


< 
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Our  method  gives  substantially  better  results  than  the  static  approxi- 
mation however. 

For  future  reference,  we  list  the  values  of  the  "constants" 

E.E.N  and  a  in  Table  2,  where 
y   V   V 

2 

E     =    <v|   ^  -  ^  +  I  -  1  +  k^|v>  /<v|v> 
^  dr         r 


2 

E     =    <y|   ^  -  ^  -  I  -  1  +  k   |y>/<y|y  > 
^  dr         r 


\  =  <yly> 


and  the  polarizibility,    a,    is 

1  /  1  I  V 

a  =  <  V  r  u  > 

N     for  k  j^  0   cannot  be   compared  directly  with  its  value  for  k  =   0»      In  order 
to  find  comparable  values  we  must   compute    \/n  /k     for  k  ^  0.      This  is 
necessary  since  the  nonaallzation  of  yv  is  proportional  to  the  normaliza- 
tion of  G  and  (5-5)    is   eqvilvalent  to 


(5-7)  7-   [sin  kr  +  tan  t]   cos  kr] 


in  the   limit  of  vanishing  k.      The   relative  normalization  of  v  and  y  is 
arbitrary  and,    in  order  to  simplify  the  numerical  work,    the  normalization 
of  V  was   chosen  such  that 


Ny=    <v|v>    =    <v^lv^>    =16.125. 
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TABLE   2. 


k  in  units  of  a 


-1 


k 

E 

y 

E 

V 

N 

V 

EX/-V 

a 

0. 

-.89 

-1.5 

.11 

.35 

l+.i| 

.2 

-.9 

-1.8 

1.9       lO"^ 

.22 

i^.5 

•  3 

-.88 

-1.9 

2.8       lo"^ 

.18 

h. 

.h 

-.8k 

-2.0 

5.^       lo"^ 

.15 

3.9 

.6 

-.72 

-2.1 

k.2        lO"^ 

.11 

3.6 

Static  VEilue  of  a  =  ^.5 


N 


16.125 


E     =       <v 


dr 


2  ^  r 


1  +  k    |v>  /<v|v> 


E„  =     <y| 


N     = 


dr 

<y|y> 


^-f  -   l  +  k^ly>/<y|y> 


1^5 


5.2  Disctisslon  of  results 

We  see  from  Figure  1  that  ouj:"  self-corLSistent  model  is  capable  of 
explaining  a  large  part  of  the  attractive  nature  of  the  hydrogen-positron 
interaction  as  found  by  Spruch  and  Rosenberg  L  J .   There  is  reason  to 
expect  that  extending  the  calculation  to  Include  spherical  distortions  would 
improve  the  phase  shifts  relative  to  Spruch' s  results.  We  base  this  on  the 
following  considerations: 

From  Table  2.,    it  is  seen  that  the  values  of  a  decrease  much  more  slow- 
ly with  increasing  energy  than  does  the  phase  shift.   This  suggests  that 
the  effective  long  tailed  potential  <xG/r  does  not  play  a  dominant  role  in 
determining  the  cross  section.   On  the  other  hand^  the  values  of 

^  \/n;  =  N/<y|y>  A 

show  a  very  precipitous  decline.   Kow^  <  v|v>  is  independent  of  energy  and 

1  I — 
consequently  the  normalization  of  y,  —  \/N  ,  is  a  measure  of  the  importance 

of  the  correlation  term.   This  indicates  that  the  strong  attraction  at  low 
energies  are  not  due  so  much  to  the  long  tall  of  the  effective  potential  as 
to  effects  of  the  correlation  at  smaller  distances.   In  Section  5  "we  argued 
that  we  should  retain  only  the  dipole  tena  since  this  makes  the  largest  con- 
tribution at  large  r;  but  our  results  indicate  that  the  most  important  part 
of  the  interaction  may  be  due  to  the  short  range  part  of  the  correlation. 
For  this  reason  we  feel  that  the  spherical  distortion  may  contribute  signi- 
ficantly even  though  it  falls  off  more  rapidly  than  the  dipole  term. 

The  experience  of  Massey  and  Moussa  L  J  lends  some  credence  to  this 
idea  that  short  range  effects  predominate.  They  introduced  a  long  tailed 
potential  into  a  variational  principle  but  they  used  a  relatively  uncorre- 
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lated  vave  function  in  comparison  to  the  one  chosen  by  Spruch  and  Rosen- 
berg L  -•  .   Despite  the  inclusion  of  polarization,  the  Massey  and  Moussa 
results  indicate  an  overall  repulsion  only  slightly  less  than  that  given 
by  the  static  model. 

The  other  ajrgument  which  indicates  the  relative  insensitivity  of  the 
phase  shift  to  the  long  tailed  potential  is  obtained  from  the  solutions  to 
(5-2,^4-)  vhen  v  is  chosen  to  be  v  .   In  this  case,  a  =  ^-5^  ^ut  the  value 
of  tan  r|/k  for  zero  energy  is  .85  as  compared  to  .92  and  a  value  of  a.  ■=  k.k 
vhen  the  full  iteration  scheme  is  used. 

It  should  be  noted  that  the  relative  success  of  v.^  in  predicting  a 

AD 

value  of  tan  r|/k  so   close  to  the  value  obtained  by  solving   (5-I)    for  v 
does  not  necessarily  mean  that  the   adiabatic  theory  vlll  give   similar  re- 
suits  since  the   form  of  y  is  given  by  QG/r     only  for  large  values  of  r. 
On  the  other  hand,    the  numerical  vork  involved  in  the   self  consistent  model 
is   considerably  simplified  with  the   choice  of  v  =  v       and,    since  the  dif- 
ference between  the  results   is  small,    it  wotild  probably  be  better  in  futiire 
applications  to  take  this   choice. 

In  this  paper  we  have   developed  a  self  consistent  Hartree  model  for  the 
scattering  of  positrons  by  atoms  and  applied  it  to  the  scattering  of  posi- 
trons by  hydrogen.      The   results  of  this   application  are  encouraging  since 
it  seems  possible  to  find  by  this  method  the  distortion  of  the   atoms  by 
the   slow  incoming  particle.      This   distortion  plays   a  jgeedomlnant  role  in 
this  problem  and  may  be  important  for  more   complicated  problems. 

While  other  techniques,    especially  the  work  of  Spruch  and  Rosenberg, 
hav«  previously  shown  the  Importance  of  short  range   correlations   in  low 
enercy  icattering,    our  method  can  readily  be   adapted  to  systematic  machine 
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calculations  of  scattering  of  positrons  by  any  atom  for  which  accurate 
unperturbed  atomic  wave  functions  are  known.   Improvements  in  accuracy 
over  the  present  work  can  be  obtained  by  Inclusion  of  monopole  and  higher 
multipole  terms. 

Application  to  the  more  interesting  and  more  important  problem  of  the 
scattering  of  electrons  by  atoms  seems  possible  although  some  of  the  diffi- 
culties due  to  exchange  have  not  been  completely  resolved.   The  simplest 
approach  would  be  to  symmetrize  our  wave  functions  after  we  have  derived 
the  coupled  differential  equations  for  the  correlation  tenns,  and  then  to 
use  the  symmetrized  wave  function  for  deteniiining  the  scattered  amplitude. 
We  feel^  however^  that  a  more  correct  approach  is  to  symmetrize  the  trial 
function  before  inserting  it  Into  the  variational  principle.   The  resulting 
equations,  however,  are  so  complicated  that  some  simplifl cations  must  be 
Introduced  before  numerical  analysis  can  begin.  We  have  not  Investigated 
in  detail  what  simplifications  are  possible.  A  resolution  of  these  diffi- 
culties would  provide  a  systematic  approach  for  the  problem  of  the  scatter- 
ing of  electrons  by  atoms. 
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Appendix  I 
Asymptotic  Expajaslons   of  G  ;,    j^  and  v 


The  behavior  of  the  integral 

In  the  neighborhood  of  r  =»  R  may  be  written  as 

R  » 

(i-i)     <g|R,lh>  =-^  /g(t)h(t)t^dt  +  /|  fi^^Mil 

o  R 


+ 
r 


i^  /.(t,.(t,t'.t  .  //  £(iMi)  «. 


If  the  last  two  terms   are  expanded  in  a  power  series  around  r  =  R,    the 
zeroth  and  first  order  terms  vanish.      Then  we  may  approximate    \g|Rylh/    as 

R  » 

(1-2)        <glR^|h>    ^  -j^   |g(t)h(t)t^dt  -.  /f    S^^  dt. 


-  "o  R         ^ 


When  R  =  0, 


(1-5)  <dB,|h>^//siiMildt. 


t 
o 


For  small  r, 

-2r 


(lA)         Vg  =  ^  (1  +  r)  ^i-  1, 

Using  (l-5)  and  (l-^)  the  equations  for  y,v,  and  G  may  be  put  in  the  form 
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(1-5) 


,2  2         *  r 

dr  r 


F(r)   =  &{r) 


F  =  V        J,    t  =  1,   upper  sign,    X  = 


F  =  y       ,    t  =  1,    lover  sign,   X  = 


F  =  G       ,    e  =  0,    lower  sign,   X 


^^I'r'^l' 


y 

+ 

2 

<y. 

^i> 

V 

+ 

2 

<'i^ 

'!> 

y 

<-i 

V 

k^ 

+ 

2. 

Solving  for  F  "by  means   of  a  power  series 


(1-6)  F  =  X]     V 


n+s 


we  have  the  reeuraion  relation. 


(1-7)     ((n  +  s)(n  +  B  -  1)   -  t(l  +   l))c  +  2c  ,  +  Xc  ^  =.  0,  n  ^  2 
^  ^     ' ^  ^     ' '   n  n-1     n-2    ' 

F  will  be  regular  at  the  origin  if  s  =  t  +  1.  With  this  value  the  first 
three  coefficients  are 


c  =  const, 
o 

c 
-     o 

1       t  +  1 


=2  = 


(2i  +  5)(t  +  1)  "  2(2i  +  5) 


One  must  take  the  inhomogeneity  into  account  in  order  to  calculate  c  . 
Clearly,  in  the  neighborhood  of  the  origin,  the  particular  integral  and 
the  complementary  function  have  the  same  form.   The  value  of  c  necessary 
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for  the  particular  integral  is  obtained  by  satisfying  the  boundary  con- 
ditions at  infinity. 

In  order  to  find  the  asymptotic  forms  of  y  and  G,  ve  note  that  u 
and  V  vanish  exponentially  for  large  r.  Then,  according  to  (l-2)  for 
large  R 


<u|r^|v>    =  ~     /     uv  t  dt. 

r      ^ 
o 


The  eisymptotlc  form  of  the  equation  defining  y  is 


R 


R 


^        2     _   2  ^  ^     +  -^  i  r  wdt  +  -^  K  Tvt^vdt 
dr^       r^       ^         y      /--,--Nr, 


■ /    wax,  -I-  — r 

<vlv>  ""J  5<vlv>  r^^ 


y 


R 

— —  ^     /      vtudtG 

<v|v>   s/5  /    { 


or  approximately, 


(1-9) 


dr^  '  r^  ""  5  <^|v>  r^ 


+  E 


2   <^l^l^>   G         ^      G_ 
^~/5     <v|v>     r^         y  r^ 


The  equation  for  G  is 


(i-io) 


dr 


G=-^    <ultlv>     ^=  7 


G     2 
r 


<v|t|u> 


7      = 


y        75         <v1y> 


^P=-4:       <u|t|v>. 
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Let 


(l-ll) 


G  =  2_J     ~H  ^^^(^^)    +   ^°^   ^  ^—d     "5 
r  r 


(1-12) 


y 


1^    ^;Xp  sin(kr)   +  cos  kr  2^    — - 


n+2 


Substituting   (l-ll)   and  (l-12)    into   (l-9)    and   (l-lO)    leads  to  the  recursive 
relations 


(1-15)  (E     -  k'')6^  +  k(n  +  l)e^_^  +   (n(n  +  l)-   2)b^_^  +  - 


<vlt2|v> 


6      ,  =  7  a 


5  <yK>     ^-5     y^ 


(1-15) 


(E     -  k   )e      -  k(n  +  l)6      ,   +   (n(n  +  l)-   2)€^   „  +  - 
^   y  '   n  ^  '    n-1        ^    ^  n-2       5 


1,  <vlt'|v> 


€      ,   =    7    P 


5     /    I    >        n-5         y'^n 


(l-l^) 


k  (3     +  n(n 
n  n 


-k  a     +  n(n 
n  n 


l)a     , 
n-1 


^G  \-5 


l)p     ,   +  7^  e      ,. 
^n-1  G     n-5 


The  leading  terms   in   (l-15)    are 


(1-15) 
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Thiis, 


_.    QG 
r 


(1-16)  Q  = 


<Cv|r|u> 


^    <v 


d         2         2       , 

dr         r 


V 


> 


From  (l-l^)  we  see  that 


(1-17)  (3^  =   p2  =  cu^  =  a^  =  0. 

In  solving  the  problem  numerically,  we  assume  the  nonnallzation  a  =  1 

o 

so  that  3  =  tan  ri. 
o       ' 

When  k  =  0,  (l-lO)  reduces  to 

(1-18)    -4  G  =  7,  ^  . 
dr        r 

The  complete  homogeneous  solution  to  (I-I9)  is  of  the  form  a  r  +  a  • 

Accordingly,  we  take 


n  T 


\r — I   a, 
n 


n+1 


which,  aTter  substitution  into  (I-I8, 11),  result  in 


(1-20)    EyS^  +  (n(n-l)-2)6^_2  +  ^  — -^  \_^   =  ^^a^ 


n(n-l)a  =  ?f  5  ^. 
^   '  n    e  n-2 
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We  take  a  =  1,  bo  that  a^    is  the  scattering  length.  From  (l-20), 

ve  see  that  none  of  the  a  vanish.  Thus  the  particular  Integral  vanlahes 

&B  —.     When  k  ^  0^  the  particular  integral  vanishes  as  — ^  .   Thus, 

r 


Particular  Integral 
Complementary  function 


,  k  =  0 


i^  ,  k  ^  0 
r 


and  ve  see  that  the  polarization  plays  a  larger  role  for  k  =  0. 
We  have  in  general  for  v 


(1-21) 


dr    r  <^y|y> 


(<y| 


Rjy>  +f<y|R2ly> 


Nov 


^<y|R^lG>u 

v5 


u  =  2r  e 


,7+1 


dt 


R  R 


^dt 
t 


r5       r2 


-  5^  - 


Thus,    for  R  large,   we   can  approximate   (l-2l)   by 


(1-22) 


•-dr 


2         ^       ^ 
-o  +  -  +  E 

2       r         V 
r 


-r      r 


<y|y>\/3  L 


R  00 

o  R     ^ 


r  A  _,  ^  2  -r    -r 
L-  +  Br     Je      . 


There   are  two  particular  Integrals  for  v, 


(1-23) 


V   = 


2-j     \     n+i  +     2-J     ^n  ~^. 


-r 


which  results  In  the  recursion  relation. 


^\  -^  ^K  +  &  +  (^  -  Dj^.i  +  [(n  -  1)  -  2]X^_ 


2     no 


(1-21^) 


(\   +  1)^^  +  0^  +  (n  -  5)]^^.^  +  [(n-M(n-3)-2]^^_2=  B6^^ 


We  keep  only  the  leading  terms, 

R 

/  yt  G  dt 
h  o 


X  = 


°    S 


<y 


d    2    2,2 


^2    2   r 
dr    r 


y> 


(1-25) 


yG 


dt 


<y 


^  _  ^  .  2  ^  j^2 

,2    2   r  ^  ^ 

dr    r 


> 
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We  can  use  this  technique  to  solve 


r  ^ 

Ldr^ 

- 

2         2 

2-^r-lJ 
r                      ^ 

1    2     -r 

by  defining 

^AD  = 

V 

-r 
e      . 

Then 

r  ^' 

Ldr^ 

- 

r 

2 

r 

V  =  -kr^   . 

Let 


•'^— J   n-£ 


which  results  In 


y  2(n-8_)_f_2 

-^-J     n-s+1     ^^—J 
n     r 


((n-s)(n.sH-l)-2)^-S;j;2=-^-' 


If  a  ^  0,    then  8=5  and.  ve  have  the  recursion  relation 


(2n  -  6)a^  +  ((n  -  k){n  -   5)  -  2)a^_^  =  - 


i^6 


no 


and 


a  =  +  1 
o 

a,  =  2a  =.2 

1     o 

a^  =  0 
2 


(I-2T)    v^  =  e- 


■'v  =  [a^r^  +  a^r^J  =  (2r^  +  r^)""'' 


e 
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It  Is  Interesting  to  note  that,    In  the   case  of  v.,^,    the  Inhomo- 

2  -r  5  -r 

geneo\is  term  goes   as  r  e        and  the  solution  goes  as   r  e      .      This  resvilts 

from  the   fact  that  the   coefficient  of  r  in  the   exponential  is  the  same 

as  the   constant  In  the  operator  and  is   comparable  to  a   resonance    effect. 

In  all  other  cases   considered  the  tvo  ntmihers   are  quite  different   and  ve 

have  the  result  that  the  leading  term  of  the  particular  integral  is   given 

by  the  inhomogenelty  mtiltlplied  by  some   constant. 

It  is  necessary  to  know  the   complementary  function  of  (l-22)  which 

is  regular  at  infinity  for  the  nttmerlcal  vork.      Thus,   we  take 


d  2  if        , 

7T  -  ^  +  r  -^  \ 
dr  r 


v(^)  .  0 


vt^'^O 


E  is  negative  for  small  k  so  we  take 


Then 


v(^)..»e-^^,    x.yis; 


,-   2 

d   .  /28 


dr 


\  ■"     '  A-r-  2       r 


2   ^r 


'dr 


z  =  0 


Let  s  =  7-  and  take  z  to  be  of  the  form 


V— 1   l3 


which  results  in  the  recursion  relation 
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2X1 


\  +  [n(n-l)   +  (f(f  -  1)    -   2)-  J(n-l)]b^_^  =  0, 


and  therefore 

2 
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Appendix  II 
Techniques  Used  In  Numerical  Integration 

Let  us  vTlte  the  equation  to  he  solved  In  the  notation 
(Il-l)    ^^^G  =   f^y 
(II-2)     o-fgV  =  f^u 

(II-3)    -^^y  =  ff 

Since  ve  do  not  know  the  boundary  conditions  on  G,  y  or  v,  it  is  necessary 
to  Iterate  the  equations  in  the  hope  that  the  solutions  will  converge. 
The  interatlon  was  started  by  using  v   for  v^   and  G^.  for  G    where 

The  functions  v^""^  and  G^°  were  used  to  solve  for  y^   ;  y    and  v^ 
for  G^''"^  G^""-^  and  y^*-"^  for  v^-^^  and  so  forth. 

The  numerical  integration  of  these  equations  was  performed  by  means 
of  a  difference  scheme  known  as  Mllnes'  Method  ^     -^  .   This  method  requires 
the  knowledge  of  four  previous  values  of  the  function  and  their  second 
derivatives  in  order  to  step  ahead.   Let  A  be  the  interval  and  suppose  we 

have  obtained  the  values  z .  ,  z'.',  1  =  1,2,5,^.   One  makes  an  estimate  of  z^ 

1   1  5 

from 

(II-J+)    z^^°^  =  2z^  -  z^  +  k^{z"   +  5^z"/5) 


where 


2_M  _ 

■'l+l  ■  "l  "^  "-l-l' 


(II-5)    6  z"  =  z_,  -  2z,  +  z 


The  second  derivative  at  1  =  5  can  be  computed  ujBlng  z    .   The  corrected 

value  of  z_  is  obtained  from 
5 

(II-6)    z^  =  -  2zj^  +  z^  +  A^(z|;  +  6^z|;/l2). 
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The  nniltlpole  Integrals 

<^|R,|v> 

vere  carried  out  by  changing  the  form  to  read 

(ii-T)       <.M.>  =  -ij  /  ^^«  .  /  /  -1^  .  /  /  -Hj  at. 

o  o  o 

The   first  two  terms   are  Integrated  from  the  origin  using  a  four  point 
difference  fornrula.      Thus 

^1+1  h 

(II-8)  f  v(t)dt  =  A    X!     CjV^_j 


x^ 


c^  =  0.57500 

Cg  =  0.79167 

c,  =  -  0.20855 
5 

cj^  =  0.0hi667 


The  Integrations   are   carried  foward  to  some  value  R  vhere  the   functions 
assume  their  asymptotic   form.      The  third  term  in  (lI-7)   may  then  be 
evaluated  by  adding  a  correction 

R     * 
R 

dt  vhich  has  already  been  computed.   Sinc6  v  and 

t  2_ 

u  fall  off  rapidly,  I   (v,v)  and  I  (v,u)  are  negligible. 

I  (yfj)    and  I  (y,  G)  are  computed  by  assuming 
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G  =»  sln(kr)  +  tan  t]    co8(kr) 

y  -  2 

r 

for  r  >  R.   I   (y,y)  and  I   (y,  G)  can  then  "be  Integrated  by  parts  a 
number  of  times  depending  on  7.   The  result  is  expressed  in  terms  of 
sine  and  cosine  integrals  the  values  of  vhlch  must  be  supplied  with 
the  initial  data.   The  definite  integrals  \v|w_)>  are  carried  out  in  a 
simileir  fashion. 

Four  values  of  the  function  are  necessary  to  initiate  the  inte- 
gration.  These  may  be  computed  from  (l-8).   The  choice  of  c  is  arbi- 
trary since  an  error  in  c  will  have  the  effect  of  introducing  the  homo- 

o 

geneouB  solution.   Thus,  suppose  we  are  integrating  the  equation 

(II-9)    ■^v  =  fz 

with  the  homogeneoxis  solution  w^        and  w^        defined  by 

(11-10)         ^J^^  =  0 


w^   =  regulax  at  the  origin 

(2) 
w^   =  regular  at  infinity. 


If  we  use  the  polynomial  (l-6)  to  start  the  integrations,  and  we  choose 

c  =  c     we  can  expect  the  trial  function  w^   to  be  of  the  form 
o    o  ^ 


(ii-n)   w^^*^  =  w  +  p.y^^  +  P^^v^^^ 


If       '^2t 

(c) 
where  w  Is  the  desired  particular  integral  with  c  =  c    .   The  coef- 
ficient, Pp„  ,  arises  from  having  truncated  (l-6)  and  can  be  shown  to 
be  negligible.   The  coefficient,  p  ,  is  proportional  to  the  difference 
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c  ^   -  c    .   It  is  also  possltle  to  start  the  integration  by  means 
of  the  asymptotic  forms  (l-12)  ,  (l-25)  and  this  vas  done  In  the  v 
Integration.   In  this  case  one  obtains 

In  this  case,  B,^  and  B_,  are  results  of  truncation  errors.   B,,  Is 
lb      da  lb 

usually  less  than  Pp,  and  since  w^   decreases  exponentially  with  de- 
creasing r,  3-1-u'w    becomes  smaller  as  the  Integration  proceeds.   (The 

(2) 
same  holds  true  of  Pp^w^   In  (ll-ll).)   In  order  to  obtain  the  partic- 
ular integral  v,  we  reqtiire  that 


s^K,J^^^v^^KsJ^K^ 


(11-12)       w^^  ^  +  p^w^  ^  \         +  h- 

Requiring  ( 11-12)  to  be  satisfied  at  two  values  of  the  argument  w,  we 
are  able  to  solve  for  p _  and  p  and  consequently  w.   An  alternative 
procedure  which  was  carried  out  for  y  and  G  is  to  carry  out  the  forward 
integration  until  the  cutoff  value  R  is  reached.  One  then  matches  w„^  ' 
with  the  asymptotic  form  (l-ll)  for  the  G  integration  and  (l-12)  for 
the  y  integration. 

In  determining  the  veilue  of  p„  for  G.^   ,  we  set  a  =«  1,  p  =0 
in  (l-n)  so  that 

(11-15)       G,,  =  sin  kr  Y"  —  +  cos  kr  3"  —  . 

AS  ^ —   n         ^ —   n 

r  r 

The  matching  equations   are  then 

(II-14)  G^^"*^^  +  pjG^''"^  =  G^  +  tan  Ti  cos  kr. 
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Reqailrlng  (ll-l^)  to  hold  at  two  points  is  sufficient  to  determine  p 
and  tan  t\. 


.(*>  1. 


The  value   of  (3     for  y  is   obtained  by  setting  (3     =  tan  r\   and 

defining  y       as   the  quotient  of   (l-12)   divided  by 


Q  = 


2 


<v|r|u> 


^     .   V 


<^ 


^2  2  ^  r       -^ 

dr         r 


■> 


We  then  require 

(11-15)  y,^'^  +  ^/^^  =  Vas- 

Again,    req\iiring  ( 11-15 )   to  hold  at  tvo  points   is  sufficient  to  determine 
Q^  and  p   .      Strictly  speaMng,    Q^  should  equal  Q  and  it  should  be  possible 
to  obtain  p„  by  matching  at   a  single  point.      In  practice,    it  is   foiind 
that  the   asymptotic  series  defining  y  does  not  give  sufficient   accuracy 
at  the  values   of  R  vhlch  were   chosen.      A  value  of  (^   slightly  different 
from  Q  will  then  give   a  better  fit  to  the   assumed  asymptotic  form. 

The  value  of  the   cutoff  radius,    R,    is   limited  by  the  fact  that 
y         — >  exp(+  \/-E     r|    and  that,    for  small  k,    E     is  negative.      If  one   is 

•J  J 

vorklng  to   an  accuracy  of  n  significant   figures,    and  exp[^\/-E     r]    is 
greater  than  10  ,    then  yj~        =   ^_^         to  n  significant  figures.      As  typi- 
cal figurea,    let  R  =  15,    Q  =  l/2,    E     =   -1 

y  ~  -^    I  •  if  X  lo"'  =  2  X  lo'^ 
R 

/^^  ^  e^5  ^  5  X  10^    . 


*'» 


tin  I 
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Thus,  y  is  contained  In  the  last  significant  figure  of  yj"   assuming 

8  significant  figures.   The  best  possible  c     must  be  in  error  in  the 

th  —ft 

9  significant  figure  and  therefore  gives  a  (3   of  the  order  of  10~  . 

The  best  p  which  occurred  was  of  the  order  of  10  so  that  y  is  re- 
duced to  6  significant  figures.  Further  reduction  in  accuracy  arises 
from  integration  errors.  Thus,  R  =  I5  is  about  as  far  out  as  one  can 
take  the  cutoff  radius. 

The  iterative  process  vas  speeded  up  by  taking  a  linear  combina- 
tion of  the  input  and  output  values  of  G.   Thus,  if  G     vas  used  to 
compute  V  and  y  and  these  functions,  v^    and  y    ,   were  used  to  com- 
pute G  ^   ,  then 
o 

.G^(1'-.X(GJ^'-G^('') 

vhere  G^      is  the  function  used  in  the  next  round  of  iterations. 
The  fastest  convergence  vas  found  for 

X  =  1.2  . 

Convergence  vas  also  aided  by  renormalizing  v  after  each  iteration. 

Strictly  speaking,  this  is  not  necessary  since  the  product  vy  is  inde- 

(c) 
pendent  of  the  normalizing  constant.   On  the  other  hand,  the  c    's 

do  depend  on  the  norma.llzation  and  vlll  be  improved  on  each  rotmd  of 

integration  if  the  normalization  of  each  function  is  kept  constant. 
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